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Abstract 

We show that the stability constant of the Lennard-Jones poten¬ 
tial in M 3 , <h(x) = |d||^ 12 — 2||x||^ 6 , is smaller than 14.316. This 
is remarkably smaller than the best previously known bound. Our 
method is very elementary, and probably applicable to other similar 
potentials such as the Morse potentials. We also improve slightly, in 
the Lennard-Jones case, the lower bound for the minimum interparti¬ 
cle distance of an optimal n-particle configuration to 0.684. 


1 Introduction 

Intermolecular forces are often well described through pair interactions 
given by a radial potential which is very repulsive at short distances, weakly 
attractive at long distances and possess a unique equilibrium distance. 
One of the most extensively considered examples of such potentials is the 
Lennard-Jones potential, [2] 4>(a;) = Ai||a ;||2 12 —A 2 ||x|| 2 6 (where A\,A 2 > 0). 
By rescaling domain and codomain we may assume, without loss of general¬ 
ity, that Ai — 1 and A 2 = 2: 

R 3 -A R $(x) = IAH 2 12 - 2||x||2 6 

which attains the global minimum at 1, and 4>(1) = —1. The term ||a;||^" 6 can 
be theoretically justified as an interaction between dipoles, while the ||x|| 2 12 
term resembles a hard-core interaction, since the potential grows fast as the 
norm of x decreases from 1. 
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For a finite configuration Q C M 3 , it is well-known that the energy per 
particle for the Lennard-Jones potential is bounded below. This condition is 
called “stability”. An introduction to this concept can be found in [5] section 
3.2, jointly with useful criteria to determine stability, e.g. proposition 3.2.8. 
An important problem associated to these kind of potentials is to estimate 
the minimum possible energy per particle. Thus, the “stability constant” is 
the minimum real B > 0 such that 

x^y 

for every finite QcK 3 . ffere \Q\ denotes the cardinality of Q. In this paper 
we prove B < 14.316 for the Lennard-Jones potential $ (theorem 5.1). 

One good reason that motivates the quest for a good estimate for B is 
its direct implication on the convergence radius of the cluster expansion of 
the corresponding grandcanonical ensemble. This expansion shows that the 
system behaves like a gas in the convergence region, and allows to com¬ 
pute exactly the thermodynamic observables. In 1963, Penrose and Ruelle 
have independently shown, for certain class of stable potentials, the conver¬ 
gence region |A| < (e 2 ^ B+1 C(/3)) _1 , where A is the activity, /3 the inverse 
temperature, and C a function that depends on $ but can be computed in¬ 
dependently from B. This classical result was recently improved by Morais, 
Procacci, Scoppola [4] and de Lima, Procacci [1], so the convergence region 
can now be written as: 

IAI < - K - 

- e pB+ ic(/3) 

where C is a function similar to C. So far, the best theoretical upper bound 
for B = Blj to their knowledge was 41.66, provided by Schachinger, Addis, 
Bomze and Schoen [6], which is significantly higher than our bound 14.316. 
Therefore, the present article enlarges drastically the proven convergence 
region for the Mayer series of the Lennard-Jones gas. 

As already mentioned in [6] , there are not many articles establishing rig¬ 
orous results on the lowest energy configurations for the Lennard-Jones or 
other similar potentials. There has been, however, a lot of effort put at com¬ 
putationally finding the optimal configurations for manageable amounts of 
particles; see for example [3], [8]. A reasonable lower bound for B can be 
obtained by considering a particular configuration: a face centered cubic lat¬ 
tice. This gives B > 8.61; see [6] for the details. Another important feature 
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is the minimal interparticle distance for the optimal n-particle configuration. 
Here we improve the lower bound from 0.67985 [6] to 0.684. A brief descrip¬ 
tion of the significance of this constant can be found in [7] section 1.1 and 
[6] section 1. 

The proof of the main theorem, the upper bound for B , follows by com¬ 
position of two main ideas. The first of them is a very simple method to 
estimate the energy of a given particle for a configuration with no particles 
closer than a given a > 0. Consider disjoint spheres of radius | around each 
particle, and define a function whose average on every ball is larger than the 
energy of the corresponding particle. Then we can just replace the sum by 
the integral over the whole space divided by the volume of each ball. 

The other idea is to estimate the global energy of a configuration with 
the method previously described, but with spheres of radius | (actually we 
need a radius a little smaller). The crucial observation is that although these 
spheres can intersect each other, we can control the contributions from the 
intersections and cancel each of them with the high energies generated by the 
two particles which are so close together as to generate such an intersection. 

In order to get good estimates, it is important to make a good choice for 
the function to be integrated. The average of a function on a ball reminds 
one of the theory of harmonic and subharmonic functions. If we have a C 2 
function with positive laplacian, then it is subharmonic, and the value at a 
point is smaller than the average on every ball centered at it. The Lennard- 
Jones potential satisfies that condition for r > ro ~ 1.14. For smaller values, 
we make a harmonic extension until a convenient radius. By the general 
nature of the arguments involved, it is reasonable to believe that the same 
strategy might be used to estimate the stability constant of other similar 
radial stable potentials such as the Morse potentials. 

Let us now summarize the structure of the exposition. In section 2 we 
introduce some auxiliary functions jointly with basic properties and formulas 
that we will need later. In section 3 we find a good bound for the energy 
of a single particle depending on the minimum distance separating any two 
particles of the configuration. As a corollary, we get our bound for the 
minimum interparticle distance for a lowest energy n-particle configuration, 
valid for every neN. We also provide a proof for the existence of minimum 
energy configurations for a fixed number of particles, a previously known 
fact. In section 4, we repeat the ideas from proposition 3.1 in order to get 
estimates involving larger balls, whose radii are almost Finally, in section 
5 we prove the main theorem, B < 14.316. 
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2 Elementary facts 

As a convention, we will work with the minus-energy, instead of the energy: 

M >0 M, h(r) = —r~ 12 + 2r~ 6 
The following function will play an important role: 

h — X(o,i] + X(i,+oo]^ 

where \R is the characteristic function of the set R , defined as 1 on R and 0 
on its complement. 

The function h has a unique root at i ~ 0.8909. It is negative for 
smaller values and positive for the rest. Its derivative h! has a unique root 
at 1, it is positive before 1 and negative thereafter. 

Formulas 2.1. Consider in M 3 two balls B\ and B 2 with radii rq and r 2 
respectively, whose centers are separated by a distance d < rq + r 2 . The 
volume of the intersection is given by 

a) | B 1 fl B 2 \ = ^(rq + r 2 - d) 2 (d 2 + 2d(r 1 + r 2 ) - 3(rq - r 2 ) 2 ) 

If we call S i and S 2 the surfaces of Bi and B 2 , the area of the portion of Si 
inside B 2 is: 

b) |Si D B- 2 1 = 7r ^"( r i + r 2 - d)(r 2 - ri + d) 

A quick way to prove these formulas is the following. Start showing that 
the height of the spherical cap SiDB 2 is equal to A(ri + r 2 — d)(r 2 — r\ + d). 
Multiplying this by 27rr 1 we obtain the area. Then we can check the volume 
formula: the derivative of the volume with respect to rq must be |Si ft B 2 \ 
and the volume at ?q — d — r 2 must be 0. 

Notation. As it was already done, for a subset X C M 3 , we denote by |A"| 
the volume, area or cardinality of A" according to the dimension of A". For a 
point x G M 3 , ||a;|| denotes the Euclidean norm. 

Observation 2.2. Calling H(x) = /r(||a;||), its Laplacian satisfies: 

A H(x) > 0 for ||a:|| > (A)!. 

A H(x) = 0 for ||a:|| = (A)l. 

A H(x) < 0 for 0 < ||a:|| < (A)l. 
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Proof. From the formula for the Laplacian in spherical coordinates, since 
there is no angular dependence, we have that the Laplacian of H at a point 
x with ||x|| = r is: 

d r (r 2 d r h(r )) = 12(—Hr -14 + 5 r -8 ) 

which is greater than 0 if and only if 5r 6 >11 and lower than 0 if and only 
if 5r 6 < 11. □ 

Proposition 2.3. 

(a) The function M >0 A M, t(r) = — ff coincides with h at 

(y)e. The first and second derivatives also coincide at that point. 

' (b) 



Proof. Part (a) is very easy to check. For part (b) we take the difference 
t(r) — h(r) and turn it into a polynomial by multiplying by r 12 . 

q(r) = r 12 (t(r ) - h(r)) = -yy'^ 12 + yyy(^y 1 )® , r 11 - 2r 6 + 1 

-L J. 1Z1 0 

We know that is a root of multiplicity at least 3. By Descartes’ rule 

of signs, the number of positive roots cannot exceed 3, so there are no more 
positive roots. Now it is easy to check that q is positive for 0 < r < (y)e 
and negative for r > (-y) 6 • □ 
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Figure: The solid line represents the function h, the dashed line represents t. 


Proposition 2.4. 


The function 




i.> 0 given by 


{ 360 H i 25 

-(— W- 

121 v 5 ; 11 

- r“ 12 + 2r“ 6 




1 

6 < r 


satisfies 

MINI) < 0(x) < 4- / f(||y||)dy 

\ n x\ Jb x 

for every i6l 3 and B x any ball centered at x such that 0 ^ B x . 


Proof. By the previous proposition, part (a), we have 9 G C 2 (M>o). Consider 
the map M^ 0 —> M>0, Q(x) = 0(||x||). Reading 0 in spherical coordinates 
we conclude that 0 € C 2 («%)• The laplacian A0 is 0 for ||x|| < (y)® and 
larger than 0 for ||x|| > (y)s, so 0 is subharmonic. 

Besides, by part (b) of the previous proposition and t(r) > 1 for r < 1, 
we have h < 9 , so: 


MIH)<9(H)<4| f 0(\\y\\)dy 

\ n x\ Jb x 


□ 
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3 Estimate for the energy of a particle and 
minimum interparticle distance 


It is useful to define, for every a > 0 

//(a) := supj^ h(x)/Q C M 3 finite, d(x,y ) > a Vx,y G Q} 

xGQ 

It is the supremum of the possible minus-energies from a point when consid¬ 
ering configurations where every interparticle distance is at least a. 

In the following proposition, we show reasonable bounds for /r(a), where 
a < 0.7. The second one is slightly sharper and its proof demands more 
effort. This improved bound is not essential for our main theorem. However, 
we included it because it allows to increase the lower bound for the minimum 
interparticle distance in an optimal configuration, from approximately 0.679 
[6] to 0.684. 

Proposition 3.1. 

I) For 0 < a < 0.7 we have: 

fi(a ) 

II) For 0.6 < a < 0.7 

y(a) 

Proof. I) Take any configuration Q such that every distance is at least a. 
For every x G Q with norm larger than 0.89, consider a ball B x of radius | 
centered at x. These balls are disjoint, therefore, applying proposition 2.4: 


24 n . . 2 , 24.05 

< — / 9(w)w dw < —r— 

a 6 J o.64 « 


24 . 9 , 26.95 

< — / 6{w)w dw < — 

^ 0.54 a 


5>w< E E 

x&Q x£Q xGQ 1 1 J±Sx 


0(\\y\\)dy< 


x£Q 
||x||>0.89 


- 171 .( 1)3 


xGQ 
cc||>0.89 

24 


9{\\y\\) d y=^i 6{r)r dr < 


26.95 


II >0.54 


'0.54 


where the integral can be solved analytically. The number 0.54 is equal to 


0.89 — 
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II) For r = ||a:|| > 0.64 + | we still can use h(x) < rj-r J B 9(\\y\\)dy but 
to prove the assertion we need, for 0.89 < ||x|| < 0.64 + | 

h(x) <7^7 [ " 064 (lly|l)4y 

\ n x\ JB x 

where 0 0 ' 64 = X(o. 64 ,+oo)$- We can integrate through spherical coordinates. 
Let c— |. For a radius w between 0.64 and r + c we must consider the sphere 
centered at 0 with radius w, and the area of its surface inside B x . By the 
formula 2.1 b), this area is equal to 

7T UJ , v , \ 7T ID , o o ox 

— (c + w — r)(c — w + r) = — {—w + Zrw + c — r ) 

Thus we have, calling A = fff(y)s 


7T 


nr+c 


25, 


0°’ M (\\y\\)dy > — / (Aw 1 - )w(—w 2 + 2rw + c 2 — r 2 )dw> 


' b t . 


' 0.64 


11 ' 


7T 

W 


•1.19 


25 


(Aw 1 - )w(—w 2 + 2 rw + c 2 — r 2 )dw 


0.64 
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We have used h(w) > t(w ) if w > (y)e (proposition 2.3) for the first in¬ 
equality, and that t(w) is positive if w < 1.49 for the second. Notice that 


1.19 = 0.89+ ^ 


We compute the relevant primitives: 


j —w 2 + 2 rw + (c 2 — r 2 ) dw = —-in 3 + rw2 + (c 2 — r2 ) w = a(w) 

12 1 

-in 3 + 2rw 2 + (c 2 — r 2 )w dw = — -in 4 + -rw 3 + -(c 2 — r 2 )w 2 = £(w) 


•1.19 


25, 


'0.64 


(Ait; 1 — — )w(—w 2 + 2rin + c 2 — r 2 )dw = 
25 

= A(a(1.19) - a(0.64)) - — (£(1.19) - £(0.64)) > 


> 0.7224c 2 - 0.7225r 2 + 1.2589r - 0.5654 


where the inequality is due to the truncation of decimal expressions. 



Restoring the factors that multiply the integral, we must prove that: 


-4- (0.7224c 2 - 0.7225 r 2 + 1.2589r - 0.5654) > h(r) (*) 

4 c 6 r 

For 0.3 < c < 0.35 and 0.89 < r < 0.64 + c. 

Claim: the left hand side is decreasing in c when 0.89 < r < 0.64 + c. 
To prove the claim, we differentiate the expression with respect to c ignoring 
the factor -jC It gives: 

—0.7224c -2 + (—3)(—0.7225r 2 + 1.2589r - 0.5654)c -4 < 0 

-0.7224c 2 < —2.1675r 2 + 3.7767r - 1.6962 

The maximum of the right hand side quadratic expression is attained at 
r = 0.8712... < 0.89. Thus, the right hand side is larger than its value at 
0.64 + c. After replacing r by 0.64 + c, we see that the remaining expression 
which is quadratic in c has the correct sign. 

0 < -1.4451c 2 + 1.0023c - 0.16692 for c e [0.3, 0.35] 

and the claim has been proven. 

Therefore, it suffices to check the inequality (X) for c = 0.35. In this case 
we have that the left hand side is 

8.343 

> 22.021 - 12.639r- 

r 

The expression is decreasing for r > 0.89. Evaluating at r = 1 we End that 
it is larger than 1, so 


o Q/1Q 

22.021 - 12.639r - --> 1 > h(r) 

r 


for 0.89 < r < 0.64 + c < 1 


Finally, 


24 2 , 24.05 

— / 0(r)r dr < — 5 — 

« J 0.64 0 




The following statement seems to be well-known. We include it for the 
sake of completeness. 
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Proposition 3.2. For every n E N there is a configuration of n points 
attaining the minimum possible energy for n points. 

Proof. The function to maximize is the minus-energy, given by M^ n A- M, 
where 

Rj 1 = {Q = (xi, ...,x n )/xi G F(Q) = ^2h(\\xi-XjW) 

#3 


The proposition follows by observing that 

sup{F(Q)} = sup{F| A '(<5)} 

for a compact K c M^ n . Indeed, we can take 

K = {Q e M^/0.65 < d(xi,Xj ) < 2 (n — 1), x\ = 0} 

To see this, we will show that any configuration has less or equal minus-energy 
than a configuration in K. 

Take a configuration Q. If there are Xi, x 3 such that d(xi, xfi) < 0.65, then 
define a as the minimum interparticle distance of Q. We have a < 0.65. Take 
x,y G Q at distance a. By previous proposition, part I), 

h(\\x — z\\) < n{a) + h(a) < —ar 12 + 2a -6 + 26.95a -3 

zeQ 

z^=x 


We can easily check that — a~ 12 + 2a -6 + 26.95a -3 < 0 for a < 0.65, since it 
is equivalent to 2a 6 + 26.95a 9 < 1, which follows evaluating at 0.65 and by 
monotonicity. Therefore, the minus-energy of x is negative, so we can improve 
the configuration by taking the particle away from the other particles. We 
can repeat the procedure until there are not two particles at distance less 
than 0.65. 

Call Q' the remaining configuration. If d(x,z ) > 2 (n — 1) for x,z G Q' 
consider the segment xz and the equispaced points x = yi, £/ 2 , •••, y n — z dr 
that segment. The distance between two consecutive points is greater than 
2. Take planes P\, ...,P n orthogonal to xz, Pk passing through y Consider 
Rk C M 3 (1 < k < n — 1) the open regions between the planes Pk and Pk+i- 
Clearly, there is at least one of those regions, say Rk, without points of Q'. We 
can reduce by 1 the distance between P k and P k +i (translating accordingly 
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the points of Q') incrementing the minus-energy provided by every pair of 
particles. Notice that through this transformation, the distance between x 
and z decreases by 1, and no interparticle distance increases. In addition, it 
does not produce a pair separated by less than 0.65 if it was not there before. 
After repeating the procedure as many times as possible, we necessarily reach 
a configuration Q" satisfying 0.65 < d(x,y ) < 2 (n — 1) Wx,y G Q". 

Finally, it is clear that we can translate so that X\ = 0. 

n 

Corollary 3.3. An optimal (lowest energy) configuration Q cl 3 of n points 
satisfies 

d(x, y) > 0.684 \/x,yeQ 

Proof. Take an optimal configuration Q, which exists by the previous propo¬ 
sition, and call the minimum distance between two particles a. Take x,y G Q 
with ||a; — y\\ = a. As seen in the previous proof, if a < 0.65, the particle 
x necessarily has positive energy, so the configuration cannot be optimal. 
Therefore, a > 0.65 and we can use proposition 3.1 II). 

0 < ^ h(\\z — x||) < /i(a) + h(a ) < —a -12 + 2a -6 + 24.05a -3 

z&Q 

Z^X 

24.05a 9 + 2a 6 > 1 

By monotonicity, we see that this does not hold for a < 0.684. □ 

4 Balls with almost 1/2 radius 

The proof of the main theorem needs to consider balls with a radius c 
as close to \ as possible, in order to avoid the factor for a small c. This 
is possible for c = 0.49, as will be shown. Then the challenge is to find a 
positive function whose average value on every ball with such a radius is 
larger than the value of h at its center (it suffices to consider ||a:|| > 0.684, 
where x is the center of the ball) and at the same time its integral on M 3 
should be as low as possible. 

q0.54 

Proposition 4.1. The function M >0 - > M> 0 given by 9 Ub = X(o. 54 ,oo)$ 

satisfies 

MINI) <fil [ ""‘“(Ibll)* 

\ n x\ JB x 
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for every x G M 3 with ||;r|| > 0.51 and B x the ball centered at x with radius 
c = 0.49. 


Proof. We divide into three regions according to ||a;|| = r G M> 0 . 

Region 1: r > 1.03. It holds by subharmonicity of Q(x) = 0(||x||) and 
h < 9, i.e. proposition 2.4. 

Region 2: 0.51 < r < 0.9. 

Again, as in the proof of proposition 3.1, we can treat the integral of 
9 °' 54 with spherical coordinates, taking advantage of the absense of angular 
dependence of 0. We use the formula 2.1 b) for the surface of a sphere inside 
another sphere. The ball B x reaches the norms 0.54 and 1, therefore: 

[ 9 0 ' 54 (\\y\\)dy > — [ (Aw~ 1 — —)w(—w 2 + 2rw + c 2 — r 2 )dw 
JB X r J o 54 11 

After a straightforward calculation as in proposition 3.1 we reach: 



25 

IT 


)w(—w 2 + 2 rw + c 2 


- r 2 )dw > —0.7558r 2 + 1.127r - 0.2516 


(the inequality is only due to truncation). Now it only remains to check 


—- —0.7558r 2 + 1.127r - 0.2516) > 1 

4(0.49) 3 r ; 


for 0.51 < r < 0.9, that can be reduced to an inequality for a quadratic 
expression. 


Region 3: 0.9 < r < 1.03 

We proceed as before, but limit the integral between 0.54 and 1.39. We 
make use of h(r) > t(r) > 0 for (y)e < r < 1.49 (proposition 2.3 (b)). 


9°- 5 \\\y\\)dy > / ' (Aw- 1 


1.39 


25, 


> —1.0199r 2 + 1.7357r - 0.5418 
and it suffices to verify 


—--^—(—1.0199r 2 + 1.7357r - 0.5418) > 1 

4(0.49) 3 r v ; 


for 0.9 < r < 1.03. 


□ 
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5 Lower bound for the average energy 


Theorem 5.1. Let M 3 A M, $(x) = ||x|| 2 12 — 2||x|| 2 6 be the Lennard-Jones 
potential. Every finite configuration QcK 3 satisfies 

t 4 E *(lk-!/ll)>-14.316 

x,y£Q 

x¥=y 

so the stability constant B is at most 14.316. 

Proof. Define 

H x (y) = h(\\y - s||), H x (y) = h(\\y - *||), Q° x 5 \y) = 9^\\\y - *||) 
for x,y e M 3 . We also consider H = H 0 , H = H 0 , @ 0 - 54 = @°- 54 . Recall 

QO.ol _ ^ 0 54 +oo )R 

We can assume that Q is an optimal configuration of \Q\ points. There¬ 
fore, by corollary 3.3 we have ||x — y\\ > 0.684 for x, y e Q. For each x G Q, 
consider B x the ball centered at x with radius 0.49. 

Let us say y ~ z whenever \\y — z || < 0.98 and y z. In this case we 
have a nontrivial intersection N yz := B y D B z . 

Take x$ G Q. 


E H »(v) < E #»(») + E( ff »w -!) 

yeQ yeQ yeQ 

y^x o y^x o y~xo 

The inequality is due to the points at a distance between 0.98 and 1 from xo- 
Applying proposition 4.1, 



y¥=xo y¥= x o 



xo^y,z&Q 


The last inequality may be intuitively clear, but a simple formal proof 
is possible. For 0 < k < \Q\ — 1, call Ek C M 3 the set of points belonging 
to exactly k balls B y with i 0 ?- y E Q. At the left, the integral on E is 
counted with multiplicity precisely k, while at the right it is counted with 
multiplicity 1 + (g), which is greater or equal than k for every integer k > 0 . 
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Now we perform the sum 


x 0 EQ 

2 Y K\\ x ~ 2/11) = E E H ^iv) < 

x,y£Q xo£Q y&Q 

x^y y^x 0 


M 

\B\ 


0O ' 54 + E(r§TE / e2- M + 2(ir(y-*j~:i)) 

2/~z I I ireQ J N yz 

x^y,z 


We will show that the parenthesis in the second term is less or equal than 0 
for every pair y ~ z. First, we can control jdy Yhx^yz In ©S ' 54 by applying 
proposition 2.4 to the balls of radius 0.684/2 centered at every x E Q. We 
have, for every w G M 3 , 


£efV) = £en*)< 

x£Q x£Q 


24 / 0 °° 9{r)r 2 dr 
0.684 3 


36 

0.684 3 


< 113 


Interchanging the sum with the integral, we find 



(30.54 

^ X 


< drll3|Af s2 


113 2 x 0.98=» (0J8 “ “WH. + 2 x °' 98 ) 
where d yz = \\y — z\\. We have used formula 2.1a) for the volume of N yz . 
With the help of a computer it can be easily visualized that 

113 1 

~Y (O8P (0J8 “ d)2(d + 1J6) + 2 ( ft W - 1) < 0 

for every d such that 0 < d < 0.98. See the appendix for the proof of this 
inequality. 

Now we can conclude 


1 

ifli 


£ Mlk-»ll)< 

x,y£Q 

x+y 


1 

2 Jb 


q°.54 _ 


12 

(0.98) 3 



< 14.316 


□ 
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Appendix. We must prove the inequality 


---—(0.98 - d) 2 (d + 1.96) + 2 (h(d) - 1) < 0 (0 < d < 0.98) 

2 (0.98) 

Define P(d) as the polynomial which results from multiplying the left hand 
side by d 12 . 

P(d) = Cl d 15 + c 2 d 13 + llld 12 + 4d 6 - 2 
113 8475 

Cl ““ 2(0.98) 3 C2 “ 49~~ 

One way would be to apply Sturm’s method to this polynomial to show that 
it has no roots between 0 and 0.98. However, we can simplify. It suffices to 
show that P' > 0 between 0 and 0.98 and evaluating at 0.98. So the problem 
reduces to 

R(d) = 15cid 9 + 13c 2 d 7 + 1332d 6 + 24 > 0 

Again, we could apply Sturm’s algorithm to this lower degree polynomial. 
But we can also efficiently find the minimum of Q in our region, since 

R'(d) = d 5 (d - p\){d- p 2 )(d - p 3 ) 

where p\ ~ —1.59958, p 2 — 0.647647, p 3 ~ 0.951934. The minimum of Q 
between 0 and 0.98 is attained at p 3 , where it is positive. 
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